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1 Introduction 

The inverse problem of the calculus of variations consists in finding conditions for the existence of 
a regular Lagrangian for a given set of second-order ordinary differential equations on a manifold, 
<t = P{QtQ), so that the given equations are equivalent to the Euler-Lagrange equations of the 
Lagrangian. In order for a Lagrangian L{q, q) to exist we must be able to find gij{q., q), so-called 
multipliers, such that 




It is shown for example in [9l [2T] that the multipliers must satisfy 



det(5fij) ^ 0, Qji = gij, 
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and conversely, if one can find functions satisfying these conditions then the equations = /* 
are derivable from a Lagrangian. These conditions are generally referred to as the Helmholtz 



conditions. The solution (gij) is the Hessian of the sought-for Lagrangian with respect to the 
velocity variables, and det{gij) 7^ is the condition for the Lagrangian to be regular. We refer 
to the recent survey [12j and the monograph [2] for comments on the history of the problem, 
milestones in the literature and accounts of the different paths that have been followed in the 
past. 

We will focus here on the case where the manifold is a Lie group. An immediate example is the 
one where the second-order system is the geodesic spray of the canonical connection on the group: 
this connection is specified in terms of left-invariant vector fields X and Y by VxY = ^[X,Y]. 
The inverse problem for this specific type of second-order system has been solved explicitly for 
almost all Lie groups up to dimension 6 by Thompson and his collaborators (see j lOl \20\ [23] 
and the references therein): in each case the authors were able to decide if a Lagrangian exists 
or not, and to provide a Lagrangian in the affirmative cases. 

The second-order equations of the canonical connection are invariant under left translations. 
Surprisingly, if a Lagrangian exists, it is not necessarily invariant. The main goal of this paper 
is to solve a type of inverse problem which is on the one hand broader than that discussed in [23] 
etc. in that it deals with any invariant system of second-order ordinary differential equations on 
a Lie group, but on the other hand more restricted in that the Lagrangian, if it exists, is required 
to be invariant also. That is to say, we will deal with the following rather natural problem: given 
an invariant second-order system on a Lie group G, when does there exist a regular Lagrangian 
for it that is also invariant under G? We call this the invariant inverse problem. The invariant 
inverse problem for the specific case of the geodesic spray of the canonical connection has been 
studied in [18]. In the current paper, by contrast, we will deal with the general invariant inverse 
problem. 

It is unfortunately not straightforward to adapt the solution of the inverse problem by the 
Helmholtz conditions to the invariant inverse problem. Clearly, if there is an invariant La- 
grangian then the corresponding multiplier matrix (its Hessian) must itself be invariant (in an 
appropriate sense). The difficulty is this: one may find a multiplier which satisfies the Helmholtz 
conditions and is invariant; one is then guaranteed that there is a Lagrangian, but not that the 
Lagrangian is invariant. Roughly speaking, to obtain the Lagrangian one has to integrate the 
multiplier, and invariance may be lost as a result. In fact extra conditions, of a cohomological 
nature, must be satisfied. The occurrence of such cohomological conditions was discussed already 
nearly twenty years ago, in a similar but more limited context, by Marmo and Morandi |13j . 
We will present a version of their result, which amounts in fact to a small generalization of it, 
in Theorem [TJ The conditions also appear in [T], using the rather different framework of the 
variational bicomplex. 

It is however possible to adopt quite a different approach from these authors, by taking advan- 
tage of invariance to carry out a reduction of the problem, which turns out to simplify it in 
concept, and to make the solution considerably more useful in applications. We next explain 
this alternative approach in a little more detail. 

Because of the invariance of our problem, G will be a symmetry group of the second-order 
system. It follows that the space of interest is effectively the Lie algebra of G rather than 
the whole manifold TG, and we can first perform a reduction. The dynamical vector field F 
corresponding to the system of differential equations, namely 

■ d ■ d 
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reduces to a vector field 7 on g given in terms of Cartesian coordinates (w*) (so tliat the are 
the components ol w & Q with respect to some chosen basis of q) by 



7 = 7' 



d 



On the other hand, if L G C°°(TG) is a regular invariant Lagrangian then its restriction to 
= TgG is a function (also called a Lagrangian) I G C°°(g). We will take optimal advantage of 
the following observation (which is proved in [14j for example, though we will give a different 
derivation below): finding a solution g(t) € G of the Euler-Lagrange equations of L is equivalent 
to finding a solution w{t) G g of the so-called Euler-Poincare equations 



(where ad* is the adjoint action of g on its dual), or equivalently if C^j are the structure constants 
of g corresponding to the basis used to define the coordinates. 



To obtain the corresponding solution g{t) of the Euler-Lagrange equations we need to solve in 
addition the equation g{t)~^g{t) = w(t). 

The invariant inverse problem on a Lie group G has therefore the following equivalent reduced 
version: if T is invariant, when does there exist a Lagrangian I G C°°{q) such that its Euler- 
Poincare equations ([1]) are equivalent to the equations = 7* for the reduced vector field 7 on 
0? As we will show in Theorem [2] below, if such a Lagrangian I exists for 7, the original vector 
field r will be the Euler-Lagrange field for some invariant Lagrangian L. The advantage of such 
an approach is that the Lagrangian being sought is simply a function of the coordinates on 
the Lie algebra g = TgG, rather than a function of the coordinates ((?*,(?*) on TG satisfying 
invariance conditions. The solution to this existence problem will be given in part by a set 
of reduced Helmholtz conditions for 7, involving a multiplier matrix (kij) which, in the end, 
is the Hessian of the function / we want to find. In addition, cohomological conditions will 
again make their appearance here. It turns out, as we will establish in Theorem [4] below, that 
one of the functions of the reduced Helmholtz conditions is to ensure that certain cochains are 
cocycles, and determine cohomology classes in the cohomology of g. What is not resolved by 
the Helmholtz conditions is whether these cocycles can be made into coboundaries; that they 
can is the additional requirement for the existence of a Lagrangian. 

The two approaches, leading respectively to Theorem [1] and Theorem UJ involve classes in the 
cohomology of g which, while differently derived, are the same. Nevertheless there is a subtle 
difference between the two forms of the inverse problem, which it is worth pointing out. The 
procedure described in Theorem [1] and [13] associates with a certain set of Lagrangians a pair 
of cohomology classes, whose vanishing is the condition for there to be a Lagrangian in the set 
which is invariant. The procedure described in Theorem U] in effect associates with a certain set 
of invariant functions a pair of cohomology classes, whose vanishing is the condition for there to 
be an invariant function in the set which is a Lagrangian. 

The geometrical framework that we will use is based on a reformulation of the Euler-Lagrange 
equations and of the Helmholtz conditions in terms of a suitable adapted frame. The requisite 
background material is given in Section 2. The solution of the invariant inverse problem using 
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invariant multipliers in the Helmholtz conditions on TG is presented in Theorem [T] in Section 3. 
The Euler-Poincare equations are derived in Section 4, and the reduced Helmholtz conditions in 
Section 5. Section 6 is devoted to the proof of Theorem [U which is the solution of the invariant 
inverse problem using the reduced Helmholtz conditions, and is the main result of the paper. 
Next, we investigate the geometric structure of the reduced space. In Section 7 we show that 
Equation ([1]) is a particular example of a so-called Lagrangian system on a Lie algebroid, where 
the Lie algebroid at hand is related in a natural way to the Lie algebra Q of the Lie group G. We 
make the link between the current set-up and Martinez's approach [T5] to Lagrangian systems 
on Lie algebroids. This will result in a coordinate- independent reformulation of the reduced 
Helmholtz conditions and of the cohomology conditions. The paper ends with some examples 
and some suggestions for future work. 

Although the paper focusses entirely on left-invariant Lagrangians, it can easily be adjusted to 
the right-invariant case. 



2 Calculus along the tangent bundle projection 

One can find in the literature several reformulations of the Helmholtz conditions that are in- 
dependent of the choice of coordinates on the manifold M: see for example O [T7]. We will 
follow closely the one given in [71 [E], which is based on a calculus of tensor fields along the 
tangent bundle projection r : TM M. By a vector field along r we mean a section of the 
pullback bundle t*TM TM, and likewise for tensor fields. A section of t*TM — > TM can be 
interpreted as a map X : TM TM with the property that t o X = t, and can be expressed 
in terms of local coordinates as 

X = X\q,q)^eXiT). 

There is a 1-1 correspondence between vector fields along r and vertical vector fields on TM. 
This correspondence is made explicit by the so-called vertical lift X^ of X, given by 

• d 
= X'^. 

Any vector field on M induces a vector field along r in an obvious way. If X = X"^ {q)d / dq^ is a 
vector field on M, its complete lift X'-' is the following vector field on TM: 

X^=X^± + ^q^±. 

dq^ dq^ d(f ' 

Here are some convenient formulae for the brackets of complete and vertical lifts: 

[a:^, y^] = [X, Yf, [A°, y^] = [x, y]^ and [a^, y^] = o. 

Here X and Y are vector fields on M throughout. 

The vertical and complete lifts require no additional machinery for their definitions. However, 
if we have a second-order differential equation field, or dynamical vector field, F at our disposal, 
say 

3g* dq^ 
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(representing the second-order equations if = /*), we can use it to define the so-cahed horizontal 
lift of a vector field along r. The horizontal lift of X G X(r) is 

Any vector field Z on TM can be decomposed into a horizontal and vertical component: Z = 
+ Y^, for X,Y E X(r). In case X is induced by a vector field on M, the three lifts are 
related as follows: 

x« = i(x°-[r,x^]). 

Another useful fact, which it is easy to establish by a coordinate calculation, is that [r,X'^] is 
always vertical. 

The Lie brackets of the dynamics T with horizontal and vertical vector fields define important 
objects for the calculus along r. It can be shown that the horizontal and vertical components 
of these brackets take the form 

[r,A:^] = + (vx)^ and [r,A:«] = (vx)« + ($(x)f . 

The operator <I> is a type (1,1) tensor field along r and is called the Jacobi endomorphism. 
The other operator, V, acts as a derivative on X(t), in the sense that for / € C°°(TM) and 
X G X(t), V{fX) = fVX + T{f)X. It is therefore called the dynamical covariant derivative. 
Finally we will need the vertical derivative associated with any X G X(r). This acts on 
vector fields along r, but is completely determined by its actions on vector fields y on M and 
on functions / on TM by the formulae D^F = and D^/ = X^(/). 

In the framework of the calculus along the tangent bundle projection the multiplier matrix is 
regarded as an operator g : X{t) x X(r) C°°{TM), that is as a type (0,2) tensor field along 
r, with local expression g = gij{q,q)dq'^ (8) dq^ . The actions of both the dynamical covariant 
derivative and the vertical derivative can easily be extended to (0,2) tensor fields along r: by 
definition, for X,Y,Z£ X(t) 

VgiX, Y) = T{g{X, Y)) - g{VX, Y) - g{X, VY) 

and 

Blg{Y, Z) = X^{g{Y, Z)) - ^(D^y, Z) - g{y, V>\Z). 

The inverse problem can now be rephrased as the search for a type (0,2) tensor field g along r 
which is non-singular and satisfies for all X, y, Z G X(t) the conditions 

g{X,Y)=g{Y,X\ 
Vg = 0, 

gmX),Y)=g{X,^Y)), 
Blg{Y,Z) = D].g{X,Z). 

These are the Helmholtz conditions in coordinate-independent form. 

It will also be desirable to have a coordinate-independent version of the Euler-Lagrangc equa- 
tions. It is easy to see that the Euler-Lagrange field F of a regular Lagrangian L is uniquely 
determined by the fact that it is a second-order differential equation field and satisfies 

F(X^(L))-X^(L) = 
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for every vector field X on M. In particular, if {Xi\ is a basis of vector fields on M then 
{Xf^X^} is an induced basis for vector fields on TM, and the following set of equations is 
equivalent to the Euler-Lagrange equations: 

T{XnL))-Xf{L) = Q. 

We now consider the effect of a diffeomorphism of M on the Euler-Lagrange equations. Let (p 
be a diffeomorphism of M and Tip the induced diffeomorphism of TM. For any X € X(M), 
r(T(^)(XC) = {TifXf and T{Tlp){X^) = {T^pXy (these are of course the counterparts of 
the bracket relations quoted earlier). Moreover, if F S X(TM) is a second-order differential 
equation field so is T{Tip)T (this is the counterpart of the fact, stated earlier, that [F, Z*^] is 
always vertical). 

Let L be a regular Lagrangian with Euler-Lagrange field F. Then T(p*L is a regular Lagrangian; 
we claim that its Euler-Lagrange field is T{Tip)~^T . The proof goes as follows. For any function 
/, vector field X and diffeomorphism 93, X{ip* f) = ip*{{T(pX){f)). We know that F is uniquely 
determined by the Euler-Lagrange equations V{X^ {L)) = X'^{L) for all X e X(M). Now 

X''{Tip*L) = Tip* i{TipXf{L)) =Tip* {r{{TipXy{L))) 

= TiTipy^r {Tp*{{TpXY{L)) = T{Tpy^T {X^{Tip*L)) . 

If L is regular and Tp*L = L then T{Tp)T = F. But although the Lagrangian uniquely 
determines the Euler-Lagrange equations, it is not in general true that the Euler-Lagrange 
equations uniquely determine the Lagrangian, so if T{Tip)T = F all we can conclude is that 
Tip*L is a Lagrangian for F; if different from L it may be called an alternative Lagrangian. 
That genuinely alternative Lagrangians (Lagrangians not differing by a total derivative) can 
exist even in the most familiar circumstances is well-known: the free particle is the most obvious 
example, and lest that look too suspiciously special we could mention also motion in a spherically 
symmetric potential in Euclidean 3-space [llj . 

3 The invariant inverse problem 

For the remainder of the paper the configuration manifold M will be a connected Lie group 
G. We will use Xg and pg to denote left and right multiplication by 5 G G. Both maps can be 
extended to actions TXg and Tpg of G on TG. 

We assume that we have a left-invariant second-order differential equation field F on TG: thus 
T{TXg)T = F for all g (z G. The question under discussion is whether F admits an invariant 
regular Lagrangian, that is, whether there is a function L on TG whose Hessian with respect 
to velocity coordinates is non-singular and which satisfies TX*L = L for all g (z G, such that 
F is the Euler-Lagrange field of L. We can conclude from the analysis at the end of the last 
section that the Euler-Lagrange field of an invariant regular Lagrangian is invariant. But if 
we start with an invariant second-order differential equation field on the other hand, and it 
admits a regular Lagrangian, then all we can conclude is that its left translates are alternative 
Lagrangians, possibly different. 

We now begin to develop the machinery we need for a deeper study of the problem. 
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By left-translating a basis {Ei} of the Lie algebra of G we obtain a left-invariant basis {Ei} of 
3C{G). Similarly, {Ei} will denote the right-invariant basis of X(G) obtained via right translation. 
These bases are related by 

Ei{g) = Ai{g)E,{g), (2) 

where {A^ (g)) is the matrix representation of ad^; in particular (e) = 6j (where e is the identity 
of G). We will identify the Lie algebra with the left-invariant vector fields: then [Ei, Ej] = C^E^ 
where the Cf^ are the structure constants of 0, and [Ei,Ej\ = —C^jEf^. (This is the convention 
in [H], for example.) 

In the following, a vector Vg in TgG will have coordinates (w^) with respect to {Ei}, so that 
Vg = w^Ei{g). Then (w^) are exactly the coordinates of the Lie algebra element w = TXg-iVg 
with respect to the basis {Ei} of g. 

The following property is true for any action of a connected Lie group on a manifold: a tensor 
field is invariant under an action if and only if its Lie derivative by every fundamental vector field 
vanishes. When the manifold is a Lie group and the action is left multiplication, the fundamental 
vector fields are the right-invariant vector fields, for which {Ei} is a basis. A function / on G is 
left-invariant if and only if Ei(f) = for all i, and a vector field X on G is left-invariant if and 
only if = 0. In particular, for the left-invariant Ej, [Ei,Ej] = 0. In view of the bracket 

relations in the two bases it follows that 

^,(4) + 4-^^ = and Af4.Gfc7 = ^™G-. (3) 

The Lagrangian L and the dynamical vector field T both live on the tangent manifold TG. To 
characterize their invariance we need to know the infinitesimal generators of the induced action 
TXg of G on TG. Given that the flow of a complete lift is tangent to the flow of the underlying 
vector field, it is easy to see that the infinitesimal generators of TXg are exactly the complete 
lifts {E^} of the infinitesimal generators of the action Xg of G on G. So a function F G G°°{TG) 
is left-invariant if and only if E^{F) = 0, and a vector field Z € X{TG) is left-invariant if and 
only if [Ef , Z] = 0. Note that Ef and Ey are invariant vector fields, by virtue of the bracket 
relations for complete and vertical lifts given earlier. The functions are also invariant; they 
are linear fibre coordinates on TG, and satisfy Ej{w^) = 6j. 

The following observations will be important. First, if a function / satisfies E]^{f) = for all i 
the / is (the pull-back to TG of) a function on G. Second, if E^ {f) = /j is a function on G for 
all i then / — /jW* is a function on G. 

Recall that we interpret the Hessian of a Lagrangian L as a type (0,2) tensor field g along the 
tangent bundle projection r : TG — > G. If T is invariant then the coefficients Kij = Ej E^ [L) = 
g{Ei,Ej) will also be invariant functions. Now when we use the Helmholtz condition approach 
to the inverse problem, if we are interested only in invariant Lagrangians we will certainly need 
to add to the Helmholtz conditions the extra condition that the multiplier g should be invariant. 
As we pointed out earlier, if we start from an invariant second-order field, it is often possible 
to find non-invariant Lagrangians. Examples of this behaviour can be found in the papers 
[10\ [20l [23] for the case of the canonical connection. The reason is that in these examples the 
extra condition about the invariance of the multiplier is usually not imposed on the problem. 
However, as we pointed out before and will shortly explain in more detail, the invariance of the 
multiplier, while necessary for the existence of an invariant Lagrangian, is not sufficient. 

The invariance of g can be defined in a coordinate-independent way as follows. We first define 
a vector field X along r to be left-invariant if its vertical lift is left-invariant. Since {Ei} is 
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a basis for X(G), it serves also as a basis for vector fields along r. Then a vector field along r, 
X = B^Ei, is invariant if Ej{E^) = 0, or if its coefficients G C°°{TG) are invariant functions. 
We will say that a type (0,2) tensor field g along r is invariant if g{X, Y) is an invariant function 
for all invariant vector fields X and Y along r. It is easy to verify that this holds if and only if 
the coefficients of g with respect to {Ei} are invariant. 

We now state and prove a theorem which shows what requirements in addition to the Helmholtz 
conditions and the invariance of the multiplier are necessary and sufficient for the existence of an 
invariant Lagrangian. Let us call a type (0,2) tensor field g along r which satisfies the Helmholtz 
conditions for an invariant second-order differential equation field T and is invariant an invariant 
multiplier for T. 

Theorem 1. An invariant multiplier for an invariant second-order differential equation field T 
determines a cohom,ology class in H^{q) and one in H'^{q). The field T is derivable from an 
invariant Lagrangian if and only if the corresponding cohomology classes vanish. 

Proof. Suppose that g is an invariant multiplier. Wc set Kij = g{Ei,Ej). By the very fact that 
we have a solution of the Helmholtz conditions wc know that there is a regular Lagrangian for 
r, say L, such that Kij = EJ E^{L). Now L need not be invariant; but from the invariance of 
the Kij we have 

= 4«(if,,) = E7^7(4°(L))=0, 

whence E^{L) = akiw^ + hk for certain functions a^i and on G. Since L is known to be a 
Lagrangian and V is invariant, 

= Ef{T{E]{L))-E'f{L))=T{a,,)-E'f{a,uw^ + bi) 
= (^Ek{aij) - Ej{aik) - au&jk) - Ejih). 

We can set to zero the coefficient of w'^ and the remaining term separately (both are functions 
on G). Prom the second we see that hi is constant. Prom the first, 

Ek{aij) - Ej{aik) - anCj,^ = 0. 

Let "i?* be the 1-forms on G dual to the Ei (so that -d = 'd^Ei is the Maurer-Cartan form, not that 
it matters); then for each i the 1-form aiji}^ is closed, from which it follows that aij = Ej{fi) 
for some functions fi on G. We have 

Ef{L) = w^Ej{fi)+hi- (4) 
note that this is of the form total derivative plus constant. Next, 

= EfE'fiL) - EfEfiL) + CtjE'={L) = w^E^ [Hfj) " E,{fi) + C^i) + C^^b^, 
from which it follows that 

Ei{fj)-Ej{fi) + Cyi = aij (5) 

is constant, and C^jbk = 0. Now we can regard the bi as the coefficients, with respect to the 
basis of Q* dual to the basis of Q with which we are working, of a linear map 6 : g ^ R, so that 
KO = biC- Similarly, the aij are the coefficients of an alternating bilinear map a : x ^ R, 
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so that a(C,r/) = aij^^^rf . We now show that, viewed from the perspective of the cohomology of 
with values in R, b and a are cocycles; that is, they satisfy the cocycle conditions 

6({e,r?})=0 and a(e, {^?, C}) + M^, {C, ^1) + MC, 1^, r/}) = 

({•, •} is the Lie algebra bracket); or in terms of the structure constants, 

hkC^j = and a^zCj^ + UjiCi, + akiC\j = 0. 

Indeed, we have just seen that b^Cf^ = 0. Operating with again on Equation ([5|) and taking 
the cyclic sum we see that aij is a cocycle too. Moreover, /, is determined only up to the 
addition of a constant; and the addition of a constant leaves h unchanged and changes q by a 
coboundary. 

If aij and hi are both cohomologous to zero, then hi = 0, and by choice of additive constants we 
can assume that Ei{fj) — Ej{fi) + C^jfi = 0. But then fi = Ei{f) for some function / on G. 

But then L — Ej{f) = L — / is invariant, and of course has T as its Euler-Lagrange field and 
has the same Hessian as L. □ 

In [13] the authors restrict their attention to Lagrangians satisfying just Ef[L) = Ej{fi), that 
is, to Lagrangians which change only by addition of a total derivative under the action of G; 
they call such Lagrangians quasi- invariant, and appeal to physics to justify this choice. From 
a purely mathematical point of view such a restriction is unnecessary, and the more general 
situation is easily analysed, as we have seen. One possible interpretation of the significance of 
the element of H^{q) is this: it is not difficult to see that hi = —Ef{£)^ where £ is the energy of 
L] so 6j = is the condition for the energy to be invariant (even though L itself might not be). 
We will have more to say about the significance of hi later. 



4 The Euler-Poincare equations 

We now turn to the reduction of T to the Lie algebra Q. 

Left-invariant functions on TG are in 1-1 correspondence with functions on the Lie algebra: on 
the one hand restriction of any function on TG to TgG determines a function on TgG = g; on the 
other hand, any function on g = TgG can be extended to a left-invariant function on the whole 
of TG by requiring it to be constant along each orbit of the action. From now on we will use the 
following convention: capital letters such as F stand for left-invariant functions, vector fields, 
etc. on TG\ the corresponding small letters such as / stand for their restrictions to TgG = g. 

A vector field Z = Ef + EJ G X{TG) is left-invariant if and only if [E^,Z] = 0, that is, 
if and only if E^CE^) = and Ef{F^) = 0. Thus Z is invariant if and only if its components 
E^ and F^ are all invariant functions. We can therefore identify them with functions and 
on the Lie algebra. Note that f^Ej\e can be identified with a vector field on TgG, since it is 
vertical; the same is not true for ^■^£'^|e, however: it is defined on T^G, but as a vector field it 
is transverse to it. 

A set {^■'} of n = dimg functions on g can be interpreted in two equivalent ways. First, the 
elements of the set could be viewed as the coefficients of a C°°(g, g)-map, namely the map 
^ : w I— > £^'^{w)Ei. A second interpretation is to view them as the components of a vector field 
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^ on 0, where ^ = ^^d/dwK This equivalence of interpretations is a manifestation of the fact 
that the vector bundles Tq — g and x g — g are isomorphic, so there is a 1-1 correspondence 
between their sections. 

The two sets {^•'} and {/-'} together define a section of the vector bundle Q x Tq — > g, or 
equivalently the bundle x g x g Q- We will adopt the following convention: an invariant 
vector field Z = E'j + EJ £ XiTG) reduces to the section z = /) of g x Tg g where 
the first element ^ = ^■'Ej is interpreted as a C°°(g,g)-map and the second / = f^d/dw^ is a 
vector field on g. In particular, for an invariant second-order field 

r = w'Ej + T^EJ G X(TG) 

the first invariance condition, Ef{w^) = 0, is trivially satisfied, so the only condition is Ef{r^) = 
0. Let A be the identity map in C°°(g,g); then F reduces to the section (A, 7) of g x Tg g, 
where 

will be often called the reduced vector field on g. 

Let L G C°°{TG) be a left-invariant regular Lagrangian with Euler-Lagrange field F. We have 
shown in Section 2 that this second-order differential equation field can be characterized by the 
equations 

F(4-(L))-4^(L)=0. (6) 

We have also shown that if L is left-invariant then so also is F. We now compute its reduced 
vector field 7 on g. 

The Euler-Lagrange equations ([6]) are of the form 

w^E^EUL) + r'ElEjiL) - Ef{L) = 0. 

With the help of ([3]), the relations between the complete and vertical lifts of elements in the two 
bases is given by 

Ef = A^^E^ + w'ClE] and ^ = A^E] . (7) 
As a consequence the first term in the Euler-Lagrange equations vanishes: 

w'E^^Ej{L) = w^AiEfEUL) = w'^EjE^iL) + w'Ai[Ef,Erm = 0. 

On the other hand, for the last term we get 

EnL) = w''ClEjiL). 

The Euler-Lagrange equations, adapted to the frame {Ei}, are therefore 

r'E^^EUL) = w'ClEjiL). 

Notice that when L is globally defined and smooth, and regular, F must vanish when = 0, 
that is, on the zero section of TG. So a necessary condition for an invariant second-order 
differential equation F to be derivable from a global regular invariant Lagrangian (or even one 
smooth and regular in a neighbourhood of the zero section) is that F should vanish on the zero 
section. In fact this is just the requirement that 6j = in Theorem [H since 

k = EnL)U=o = F^.=o(A"(^)). (8) 
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Let I € C°°{q) be the restriction of the left-invariant Lagrangian L G C°°{TG) to the Lie algebra. 
Then the restriction of E^{L) to g is dl/dw^ , and so on. The defining relation for the reduced 
vector field 7 G of L is therefore 

These are the so-called Euler-Poincare equations jH]. 

Evidently if I is globally defined, smooth and regular on Q then 7 must vanish at the origin 
(this is the counterpart of the property of T noted above). So for a vector field 7 on g to be 
derivable via the Euler-Poincare equations from a smooth and regular (reduced) Lagrangian it 
is necessary that 7(0) = 0. We will come back to this point later. 

The Euler-Poincare equations should be interpreted as differential equations with solution w{t) 
in the Lie algebra. We have chosen the coordinates (w*) in such a way that they are not only the 
coordinates for w = w'^Ei in g, but also the fibre coordinates of any translate Vg = TXgW G TgG. 
To find the solution {g{t),g{t)) G TG of the Euler-Lagrange equations that corresponds to w{t), 
one simply needs to integrate the equation g~^{t)g{t) = w{t). 

So far as the inverse problem is concerned, we can use the foregoing analysis to reduce the 
problem to one on g, as set out in the following theorem. 

Theorem 2. Let T he an invariant second-order differential equation field on a Lie group G, 
and 7 the corresponding reduced vector field on Q . Then T admits a regular invariant Lagrangian 
L on TG if and only if ^ admits a regular Lagrangian on Q, in the sense that there is a smooth 
function I whose Hessian is non-singular, such that 7 is the vector field uniquely determined by 
the Euler-Poincare equations of I. 

Proof. Clearly, if L is a regular invariant Lagrangian for T, its restriction / to g is a regular 
Lagrangian for 7. Conversely, suppose that / is a regular Lagrangian for 7 on g, and let L be 
the unique invariant function on TG which agrees with / on T^G = g. Consider the functions 

^,=r'^E]^EnL)-w'ci,Ej{L), 

where F = vu^E^ -\- T^E^. We showed earlier that is invariant, and so is vu^. Since E^ 
commutes with EJ, both EJ(L) and E^EJ{L) are invariant. So ipi is invariant. But the 
restriction of ifi to g vanishes, by the Euler-Poincare equations; so ipi vanishes everywhere on 
TG. But as we showed earlier, the vanishing of (pi is equivalent to the Euler-Lagrange equations 
for L. Moreover, L is regular since / is. Thus L is a regular invariant Lagrangian and T is its 
Euler-Lagrange field. □ 

5 The reduced Helmholtz conditions 

In this section we will show that in the case of an invariant Lagrangian, not only the Euler- 
Lagrange equations, but also the Helmholtz conditions can be restated as conditions at the level 
of the Lie algebra. 

Recall that we interpret the Hessian of a Lagrangian as a type (0,2) tensor field g along the 
tangent bundle projection r : TG G. Due to the invariance of the Lagrangian the coefficients 
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Kij = EyEj(L) = g{Ei,Ej) will also be invariant functions. In what follows we will denote the 
restrictions of these functions to g by kij . 

Let us now evaluate the Helmholtz conditions, which we have stated in a coordinate free way in 
the second section, in the basis {Ei}. The first conditions are simply 

detiKij) ^ 0, Kij = Kji. (10) 

The Jacobi endomorphism and the dynamical derivative are determined by the horizontal struc- 
ture on TG. Since [r,^''] = -E^ + {w^C^i - EJ{Y''))El, it is easy to see that the horizontal 
lift of Ei is 

Ef = Ef + \ {-w^C^i + ir(r')) El = Ef - K^iEl 
say. Now both Ef{w^) = and EfiJ^^) = 0, so all Ef are invariant. From now on 

denotes the restriction of the invariant function to g. It is easy to see that the horizontal lift 
of an invariant vector field along r is invariant, and vice versa. 

We next consider the dynamical covariant derivative V. We have [F, EJ] = —Ef + {S/EiY. Now 
both r and EJ are invariant, so by the Jacobi identity [F, E^] must be invariant also. Since the 
horizontal part of the bracket, —Ef, is invariant, (V-Ej)^ and therefore {VEi) must be invariant 
in turn. In general, if X = X'^Ei G X(r) is invariant, then VX = X^VEi + F(X*)£'j is also 
invariant. We may summarize this result by saying that V itself is invariant. Furthermore, the 
coefficients of V with respect to the invariant basis are invariant functions, which can be reduced 
to functions on Q. In fact we can calculate [F, £^>'] explicitly, obtaining [F, £'>'] = —Ef + A^E^, 
so that 

VEi = iKCj=, - ^/(F'^))^fc = A^Ek, 

and the coefficients are just the functions which we know already to be invariant. 

Given that E^{Kij) = the Helmholtz condition = 0, when evaluated on the pair {Ei,Ej), 
gives 

r'^EliKij) - KkjA^ - K,k^ = 0. (11) 

The components of the Jacobi endomorphism with respect to the current basis can be calculated 
from [F,£^?]. One finds that 

^{Ej) = {^rEjE]{v') + \rc\j-\E]{r)EU^') 

- Icy^EliV') + \w'^C\kE]{T'') - \w'''w'^Ci^Cit) Ei = 
Again, the coefficients are invariant functions, and restrict to functions on g given by 

= 27 Q-iQ^ + 27 - - -,CijW ^ + 4^ - 4^ ^ C^fink- 

This somewhat uncouth-looking formula can be civilized by expressing it in terms of the quan- 
tities 
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when it becomes 

Again, for any invariant X, ^{X) is an invariant vector field along r. The Helmholtz condition 
involving the Jacobi endomorphism is simply 

Kij^i = Kik^. (12) 

Finally, the D^-condition is 

Er{K,,) = EnKij). (13) 

The conditions (jlOp . (jlip . (jl2p and (jl3p are all invariant; it is therefore enough to find a solution 
/cjj € C°°(g) of the restriction of these conditions to = TgG, which may be called the reduced 
Helmholtz conditions. The solution of the full conditions on TG can then be found by left 
translating the solution on g. 

For any 7 = j^d/dw^ £ ^(fl); we call a matrix (kij) of functions on g a multiplier matrix for 7 
if it satisfies the reduced Helmholtz conditions 

(iet{kij) ^ 0, kij = kji, 
dkjj _ dkij 

We have shown 

Theorem 3. Suppose given an invariant second-order differential equation field T, with reduced 
vector field 7. Then there is an invariant multiplier matrix (Kij) for T on TG if and only there 
is a multiplier matrix (kij) for 7 on g. 



6 The reduced inverse problem 

Theorem [2] shows that the problem of finding an invariant regular Lagrangian for an invariant 
second-order differential equation field on TG can be reduced to that of finding a regular La- 
grangian for the reduced vector field on g. From Theorem [3] we can infer that the existence 
of a multiplier matrix, that is, a solution of the reduced Helmholtz conditions, for the reduced 
vector field on g is a necessary condition for it to admit a Lagrangian. However, as we know, the 
relationship between Helmholtz conditions and Lagrangians in the invariant inverse problem is a 
little more complicated than is the case for the ordinary inverse problem. While the existence of 
a multiplier matrix on g is sufficient to guarantee the existence of an invariant multiplier matrix 
on TG, the existence of an invariant multiplier matrix on TG is not sufficient to guarantee the 
existence of an invariant Lagrangian on TG. The following theorem supplies in effect the extra 
conditions, working now entirely in terms of reduced quantities on g. 

Theorem 4. A multiplier matrix for 7 G X(g) determines a cohomology class in H^{q) and 
one in H^{q). The vector field 7 is derivable from a Lagrangian if and only if the corresponding 
cohomology classes vanish. 



13 



Proof. Suppose the functions kij on g satisfy the reduced Helmholtz conditions, so that (kij) is 
a multipUer matrix. Prom the last of the Helmholtz conditions, 

dkik dkij 

and the assumed symmetry of kij in its indices, it follows that there is a function I on such 
that 

dw'^dw^ ' 

I is determined up to the addition of a term linear in the w'' (and the addition of a constant, 
but this we can ignore). Then 



d 



dl 



Cijw'^ku. 



Let us denote the term in brackets on the left-hand side (whose vanishing is the Euler-Poincare 
equations) by Vj. Then the Helmholtz condition 



is equivalent to 



It follows that 



+ 



J _ 



dw^ dw 



0. 



whence 



= 0. 



dw^dw'^ 

But this says that there are constants /lij and Ui, the /lij being skew in their indices, such that 



(14) 



As before, we can regard the Vi as the coefficients of a linear map : g ^ R, and the jiij as the 
coefficients of an alternating bilinear map : g x g ^ R. We now show that u and /x satisfy 
the cocycle conditions t'/C]^- = and fiuCji^ + fijiCl- + HkiC\j = 0. In fact these conditions hold 
by virtue of the Helmholtz conditions, especially the condition ku^^j = kji<j)\. 

Now jjiij is half of the skew part of 



udk 



so that 



Earlier, we set 



1 fd-/ 



+ C% w'-^ ] hi - Cf 



dl 
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we now put 

II k 
Xij = i^ikji — i^jkii = + . 

We now consider (p^kji — (f^^jku, where (as we showed ear her) 

We look first at the terms in (pikji — (p^jkii which involve 7('i/'i): these are 

7(V'i)fcj« - 7(^j)^iZ = liXij) - 4l{kji) + i'jlikii). 

We substitute for the "){kii) terms from the appropriate Helmholtz condition, and find in the 
end that 

(t>\kji - (t>\ku = -fixij) + w^CiiXji - w^CijXii = 0. 
Now since fiij is constant, 



so that 



= l{X^J) + {X^A + X]lCi,)w'' 



I I I I k I I I k 



Since this expression is affine in with constant coefficients, these coefficients must vanish. 
Therefore both C\jVi = and fJ-uCj)^ + fJ-jiC^j^ + fJ-kiClj = 0, as required. 

If we change / to /' = / + 9kW^, the corresponding change in the cocycles is from (z^, /x) to 
where = i^i and /i-^- = /Ujj - OkCip oi v' = v and fJ.'iCv) = KC,v) - (^i{C,v})- That 
is, both components of (z^, u') and (/u, fi') belong to the same cohomology class, respectively. 
If the cohomology classes of u and fi vanish then we can find 9 such that /' = I + 9kW^ is a 
Lagrangian. □ 



By setting = in Equation ()14p we see that 

But as we pointed out earlier, it is a necessary condition for 7 to be derivable from a Lagrangian 
/ on that 7(0) = 0. The significance of the vanishing of z^ as a condition for 7 to be derivable 
from a Lagrangian is clear. 

We have derived two sets of conditions for the existence of an invariant Lagrangian, each involv- 
ing a pair of cohomology classes. One would hope that the two pairs of cohomology classes are 
the same. This is in fact the case, as we now show. 
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First we show that bi and are the same constants. From Equation ([8]) we have 

Since bi is constant it is enough to evaluate the right-hand side at e; here the distinction between 
Ejy and disappears, and we obtain 

To find the relationship between aij and /ijj it turns out to be convenient to work entirely in 
terms of the right-invariant fields Ei] in the end we will evaluate everything at e, using the 
constancy of the ajj, and again we can take advantage of the fact that at e the distinction 
between Ei and Ei disappears. The relations between the complete and vertical lifts of Ei and 
El given in Equation ([7]) now come into play. In particular, F = AlEf + T^AjE^; we set 
A)w^ = v\ A'p = P. 

It follows from Equation dH) that Ei{fj) = EJ{Ef{L)), whence 

E,{f,) = EnEfiL)) = Enr{Ej{L))). 

Using the expression above for F, and the evident fact that EJ{v^) = J^, we find that 

EnEj{L))+v'EnE^AEj{m 

+ Enf')El{E]{L)) + T'EliEjiEjiL))) - Ej{E'^{L)) = 0, 

or (remembering that [Ei,Ej\ = —CijEk) 

Er{Ef{L)) - E]{Ef{L)) + C^,El{L) 

= TiEUEjiL)) - Clv^Ej{E]{L)) + EU^'WkiEj (L)). 

On taking the skew part we find that 

EriEfiL)) - E]{Ef{L)) + C^^EliL) 

= \ [EUt') + Ctv') El{E]{L)) - \ {e]{V^) + d^v') EUEni)). 

The left-hand side is aij + Cfj{E^{L) — fk). At e, the right-hand side is 

1 f9^^ , ^fc„J 



2 

Thus 



kjk - I + Cj'jw^'^ hk = ipikjk - ''Pjkik = Xij- 



a,i = Xi,-C':^{El{L)-h)\,. 

Now let / be the restriction of L to T(,G. Of course L is not assumed to be invariant, so this 
differs from the association between / and L given earlier; nevertheless, it is true that 

where [kij] satisfies the reduced Helmholtz conditions. So we can write Uij = fiij + C^jfk{e). It 
is apparent that aij and fiij define the same cohomology class (they differ by a coboundary). 
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7 The Lie algebroid 



Our policy while working on TG in earlier sections was to write everything in terms of G- 
invariant quantities, that is, quantities determined by their values on T(,G = Q. This paves the 
way towards expressing the whole theory in terms of 0, or more accurately in terms of a vector 
bundle over g, namely g x Tg — > g. We can identify invariant vector fields on TG, via their 
restrictions to TeG, with sections of g x Tg g, as we pointed out earlier. The bracket of two 
invariant vector fields remains invariant, and so the bracket of vector fields on TG determines 
a bracket of sections of g x Tg — > g. This is evidently R-bilinear and skew, and it satisfies the 
Jacobi identity by construction. We will now obtain an explicit formula for this bracket, and 
deduce that it is a Lie algebroid bracket, i.e. a bracket with the above properties that satisfies an 
appropriate Leibniz rule when sections are being multiplied with functions on the base manifold. 

Let i^Ef + X^EJ and rfE"^ + Y^EJ be two invariant vector fields, so that ^{ff) = Ef(X') = 
E'j{rf) = E'jiY'^) = 0. These invariance conditions, when expressed in terms of the vector fields 
of the invariant basis, become for example (^*) = w'^CljE^ {^'^), using Equation ([7]). Thus 

[eEf,rfEf] = (rVG^- + Cw^C^jiEn^') - r^' w'^ C],Ej [e]) E^ 

while 

[eEf,Y^Ej] = -Y^E]{e)H + [ey'cf^+ew^c]M{y'')) ei. 

The bracket may be written as follows. We identify ^, r/ with g-valued functions on g, X, Y with 
vector fields on g; ^ is the vector field corresponding to ^. We think of w^C^j as the components 
of a type (1,1) tensor field on g which we denote by A: thus 

d 

The Lie algebra bracket {•, •} extends naturally to an algebraic bracket on g-valued functions 
on g, so that {^,r/} = S^^rj^Gj^Ei. Then 

[(e,x),(r/,y)i = {{^,^} + A{0{v)-Am)+x{v)-Y{0, 

[A{0,Y] - [A{n),X] + AiYif}) - A(X^) + [X, Y]) . 

For any function / on g we have liC,X), f{rj,Y)j = m,X),{r],Y)j + p{^,X){f){^,Y), as 
required, where the so-called anchor of the Lie algebroid is given by 

p{^,X)=A{0 + XeX{e). 

Thus the bracket [•, ■] does indeed define a Lie algebroid structure on g x Tg — > g. 

We denote by the section {Ei, 0) of g x Tg — g, and Wi the section (0, d/dw^); then {cj, Wi} 
is a basis of sections, and we have 

We denote by b the induced exterior derivative operator on sections of exterior powers of the 
dual of the algebroid, and by {e*, VF*} the basis dual to {e^, Wj}. Then for any function / on g. 
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while 



Using these formulae we can express the Euler-Poincare equations in terms of the Lie algebroid 
structure, as follows. The vertical endomorphism S on the Lie algebroid is just X) = (0, ^). 
The invariant Lagrangian is represented by a function / on g. We define, in analogy to the usual 
case, a Cartan form 9 and an energy function £ by 

e = S{6l), £={A,Sl)-l, 

where A = w^Ci and A = w^Wi. As we will show by a direct calculation, provided that I is 
regular (in the sense that its Hessian is non-singular) the equation 

i^se = -6£ 

determines a unique section T, which is of second-order differential equation type, so that it 
takes the form w^Ci + 7*Wj, and the 7* satisfy the Euler-Poincare equations for I. In fact 

91 i 



dw 

69 = 



( -^T^^^'Cl - i#rC^l A - ^r^^e' A 



£ = w'—-l 

Let us write T = ^*ej -|- PWi. Then the vanishing of the Wi component of irS9 + 6£ gives 

( C + ) QyjlQy^j 

whence = w'^ when / is regular. When this result is inserted in i-p69 + 6£ several terms cancel, 
and the remaining terms in the ej component reduce to 

j d I ^ dl 

as required. 



The above derivation of the Euler-Poincare equations was inspired by Martinez' framework [15] 
for Lagrangian systems on a Lie algebroid. This framework is based on the so-called prolongation 
algebroid of the underlying Lie algebroid. It should be remarked that although the underlying 
algebroid of the current system is just the Lie algebra 0, the algebroid we have defined in this 
section does not coincide with the prolongation algebroid of the Lie algebra. The prolongation 
algebroid can most easily be defined as follows. Observe that both components of a section of 
X Tq — > have a natural bracket structure. For the first, we have the natural extension of 
the Lie algebra bracket to C°°(0, 0)-functions, and for the second component we have the Lie 
bracket of vector fields. The easiest way to combine these two into one Lie bracket structure is 
as follows: 

X), (r?, Y)f = ({e, r/} + X{rj) - Y{0, [X, Y]). (15) 
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It is easy to check that this is a Lie algebroid whose anchor map p : Qx Tq Tg is simply the 
projection on the second component. In our basis 

lei, ejf = C^jCk, lei, Wjf = 0, [W^, Wjf = 0. 

A short calculation reveals that indeed the expression irS^O^ = —6^£ leads again to the Euler- 
Poincare equations. It is, however, easy to guess the relationship between the two algebroid 
structures: the section map {^,X) i— > {^,A{^) + X) is an isomorphism of the first of the Lie 
algebroids with the second. The more complicated structure of the first Lie algebroid (or of that 
presentation of the common Lie algebroid if one regards isomorphic algebroids as identical in 
principle) arises from our desire to work always with invariant objects on TG and objects on Q 
derived from them by restriction. 

We conclude this section by showing that the restrictions of the horizontal lift, the Jacobi 
endomorphism and the dynamical derivative to Q have a direct interpretation in the first Lie 
algebroid. 

At the level of the Lie algebra, 0- valued functions on Q (sections of g x — > play the role that 
the vector fields along the tangent bundle projection played at the level of TG. That is, there 
is a well-defined vertical lift of such a function = ^'^{'w)Ei to the section = ^^Wi of the Lie 
algebroid. Moreover, one can easily verify that each vector field 7 on g defines a horizontal lift 

e« = - A^W^) 

to sections of the Lie algebroid, or equivalently, a splitting of the short exact sequence 

O^{O}xr5^0xr0 ^5X0^0 

of vector bundles over g. Observe that the restriction to g of the horizontal lift of an invariant 
vector field along r is in fact the horizontal lift of the restriction to g of that vector field along r. 

The functions we have introduced in previous sections have a nice interpretation in the 
algebroid set-up. The projection by the anchor map of the horizontal algebroid section Ef = 
Cj — \lWj to a vector field on g is exactly 

p(Ef) = (^^C7L-Ai)^ = ^^A 

As before, we can define a Jacobi endomorphism and a dynamical derivative by considering the 
horizontal and vertical parts of the brackets of the algebroid section F = tu^Cj + 7*Wj: 

ir, r?l = -rf + iVr,r, [L, r?«] = (Vr?)« + ($(17))^, 

where V acts like a derivative in the sense that for / G C°°(g), V(/7y) = /Vr? -|- j{f)r], and $ 
is tensorial with coefficients as before. We have VEi = X^E^, and in fact in both cases the 
operators are nothing but the original operators restricted to g = T^G. 

We define the vertical derivative for ^ € C°°{q,q) as the map : C°°(g,g) C°°(g,g), 
determined by D^/ = ^(/) for / G C°°(g), D^?? = for a vector 77 of g (in other words a 
constant element of C°^(g,g)), and the obvious Leibniz rule for multiplication by functions. 
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We can now restate the reduced Helmholtz conditions in a coordinate-free form. The multiplier 
(kij) is a matrix of functions on g, or equivalently a map k : C°°{Q,g) x C°°{q,q) — > C°°{q), 
which satisfies the conditions 

det/c/0, ki^,7])=kir],0, 
Vk = 0, 

kmv),o = k{v,HO), 

Bjk{rj,0 = ^;k{^,0- 

for alU,r?,C e C°°(S,S). 

Suppose that the Hessian of / G C°°{g) is k. Then if ^ = {dl/dw^)e^ as before, we can define a 
2-form 6^6 on the algebroid by requiring it to vanish whenever one of its arguments is a vertical 
section (i.e. it is semi-basic) and by setting 

Then 6^9 = ^fJ-ije^ A , where 

/X,,- = p{Ef){e{Ef))-p{EfmEf))-9{lEf,Ef}) 

dl 

= ^kij-^'jku-C^^-^. 

These coefficients are exactly those we have encountered in the previous section. Recall from 
the proof of Theorem 4 that the reduced Helmholtz conditions ensure that the p,ij are constants 
and that they form a cocycle. A necessary condition for a Lagrangian to exist is that /Xij is a 
coboundary. We can now re-express this statement in terms of the Lie algebroid. For a 2-form 
fi = 5^6 with constant coefficients we get that 5fi = —^fHjCi^e'' f\e^ f\eK Therefore, (5-closure of 
the 2-form /i amounts to the cocycle condition. On the other hand, the condition that the ^i^j are 
of the form otkC^j for some ctfc is equivalent to jj, being exact. Similarly, the reduced Helmholtz 
conditions ensure that the semi-basic 1-form v = iySO + 5£ — irp. has constant coefficients 
that form a cocycle, or that 6iy = 0. Theorem 4 states that v should vanish for a Lagrangian to 
exist. 



8 Examples and applications 

The method of reduced Helmholtz conditions really comes into its own when one has to deal with 
any specific problem. In practice — certainly, if the following examples are representative — the 
cohomological conditions do not play much of a role. Where there is no invariant Lagrangian 
this is because the reduced Helmholtz conditions fail, often at the level of regularity. Where one 
is able to find a solution of the Helmholtz conditions one is usually able to integrate it by hand, 
and check directly for which integration constants the Euler-Poincare equations are equivalent 
to the equations associated with the vector field 7. 

To save space, in the following examples we will write kiji for dkij/dw\ and we will implicitly 
assume that the conditions kij = kji, kiji = kuj and so on are satisfied. 
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8.1 The canonical connection on a Lie group 



The canonical connection on a Lie group is defined as a covariant derivative operator by VxY = 
^[X,Y], where X and Y are any two left-invariant vector fields on G. As we mentioned in the 
Introduction, the invariant inverse problem for the canonical connection has been studied by 
Muzsnay in [18]; however, he uses methods different from ours. 

The connection coefficients of the canonical connection with respect to the left-invariant basis 
{Ei} of X{G) are just ^Cjk- So the coefficients F* of the corresponding second-order differential 
equation field (the geodesic spray) are in this case = ^C^i^w^w^ = 0. The reduced equations 
are therefore simply = 0. In fact the geodesies through the identity of G are just the 
1-parameter subgroups. 

If a left-invariant Lagrangian L exists, then 

r)l 

C^.4-(LK = or G^^^w^ = 0. (16) 

In view of relation d?]), L must also be right-invariant, i.e. E'j{L) = 0, and thus bi-invariant. At 
the level of the Lie algebra, this means that / G C°°(0) will be ad-invariant. This observation is 
in fact Proposition 2 in [18]. Thus a Lagrangian is a function which is constant on the adjoint 
orbits in Q and whose Hessian is non-singular. 

We will use our methods to investigate the invariant inverse problem for the canonical connection. 
The first observation is that in this case the reduced Helmholtz condition kiicl)^^ = kji(f)[ is a 
consequence of the other conditions. Since 7* = 0, 

8k- ■ 

- fefc.Af - k.kX'; = \w\kk,Ct + hkCf^) = 0. (17) 
On the other hand, 0'. = -jW^w^-G^jGIj^. But 

wVGi^^GUku - wVCifilkk.i = -wVCifilM + wVCi,Gl^kki = 0, 

so the condition kii(p^j = kji4>\ holds by virtue of condition (jl7p and the symmetry of kj^i. 

A second general remark concerns the cohomological conditions. In this case Equation (jl4p 
reads 

1 k 

from which immediately I'i = 0. Moreover 

81 

so is a coboundary. Thus the cohomological conditions are automatically satisfied for the 
canonical connection. Any function / whose Hessian satisfies the reduced Helmholtz conditions 
and is such that G^jdl/8w^{Q) = will be a Lagrangian; in particular, if I satisfies the reduced 
Helmholtz conditions and we set 

then /' will be a Lagrangian. So the inverse problem for the canonical connection reduces 
essentially to the analysis of condition (fTTl) . in the form w\kkjCf^ + kikCfj) = 0, and the condition 
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f^ijk = kikj- Where there is no Lagrangian this will often become apparent by the fact that there 
is no non-singular (kij) satisfying the first of these conditions. 

We will examine two specific situations, one in which there is no Lagrangian, one in which there 



IS one. 



The first case is that of the Heisenberg algebra, which is a 3-dimensional algebra with the single 
non-trivial bracket relation {EijE-^} = E2. Condition ()17p amounts simply to 



Evidently ki2 = ^22 = ^32 
Helmholtz conditions. 



ki2Uf^ = k22W^ = k22W^ = k32W^ = -kl2W^ + /i;32W^ = 0. 

0, and there is no non-singular 3x3 matrix 



satisfying the 



For our second example we take the 4-dimensional Lie algebra with bracket relations 

{E2, E'i} = El, {E2, £"4} = E2, {i?3, E4} = —E3 

(this is the algebra A4^8 in the classification of Patera et al. [I9j). Condition ()17p says in this 
case that the matrix 
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k22 


k23 


k2A 







4 

w 





—w 


ki3 


k23 


k33 












4 

—w 




ku 


k2i 


ksi 


^44 

















must be skew-symmetric. This leads to the following 7 independent equations for the 10 un- 
knowns kij (with i < j): 



kiiw + ki2W 
ki2W^ + k22w'^ 
k24w'^ + ks^w^ 

k22'Up' - k2Aw'^ 

k33W^ - ksiw"^ 



= kiiw'^ + ki3W , 
= ki3w'^ + k33W^, 

0, 

(ku + k23)w^ 
(ku + k23)w^- 



Evidently /C44 is unconstrained by these equations. It turns out that A;24 = ^34 = 0. The 
remaining unknowns can conveniently be expressed in terms of ku and A;23. If for convenience 
we set fell = {w^)'^F (for w'^ / 0), ^23 = G and ^44 = H then {kij) is 



-w^w^F G 



-w'^w'^F ViPur'F — G 



w'^w^F - G 
We next look at the conditions kijk 

From A;224 = ^242 = we obtain 







G 

{w^fF 




kikj. From A:i24 = ^241 
,d{w^F) 





H 



we find that 



w 



3\2 



dF 



0. 
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It follows that F = 0, except possibly where w'^ = ov w'^ = 0. Thus ku = 0, except possibly 
where w'^ = or w'^ = 0; but then by continuity ku = everywhere; and similarly for the other 
coefficients involving F. We are left with 












-G 








G 








G 








-G 








H 



This is evidently non-singular provided that G is non-zero, whatever H may be. Continuing to 
analyse the consequences of the condition kij^. = kikj we find that G must be constant and H 
must be a function of alone. This gives as potential Lagrangians 

where A and the as are constants with A non-zero, and h is an arbitrary smooth function of 
its argument. According to the general remarks made earlier, I will in fact be a Lagrangian if 
and only if qi = Q2 = as = {h doesn't play a role here because Cfj = 0). It is easy to check 
this directly. In fact for the potential Lagrangian above it is easy to see by direct calculation 
that all Ui and almost all ^ij vanish, except that /i23 = — ai, /i24 = —02 and /U34 = 03 (and 
their skew counterparts). We have shown that there will exist a Lagrangian V = 1 + OkW^ whose 
Euler-Poincare equations are exactly the equations associated to 7 if we can find 9k such that 
A^ij = GkC^y One easily verifies that this condition is only satisfied for 0^ = — a^. The sought-for 
Lagrangian /' is therefore the one above where one sets a^ = 0. 

It is interesting to note that the most general Lagrangian in this case is not just a quadratic 
form. 



The method used by Muzsnay in [18] deals directly with the equation 

as a set of partial differential equations for /. In effect, Muzsnay derives an integrability condition 
for this equation by differentiating it, to obtain 

The part of this equation symmetric in i and j is our Helmholtz condition (jl7p . the skew part 
states that the cocycle /ijj must vanish. The examples we have considered above are two of the 
many examples dealt with in p!8]. Muzsnay's results are of course broadly the same as ours; 
however, in the second case though he shows that a Lagrangian exists he does not indicate 
how to find one, whereas we have obtained the most general one. As Muzsnay points out, the 
example is also treated in [10]. By using only the unreduced Helmholtz conditions on TG, the 
authors of [10] look for a (not necessarily invariant) Lagrangian L for the canonical geodesic 
flow on any 4 dimensional Lie group. Although they are not able to give an expression of the 
most general Lagrangian, they observe in the case of the Lie algebra A4^8 that the quadratic 
part of the Lagrangian above (written in terms of invariant forms on G in their set-up) generates 
the flow of the canonical connection. They also notice that the quadratic part is a bi-invariant 
metric (as the theory predicts). 
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8.2 The Bloch-Iserles equations 



These equations appear in e.g. [3] l^- The space of interest is Sym(n), the Unear space of 
symmetric n x n matrices. The equation is 



w 



[w^N], 



(18) 



where w € Sym(n), is a skew-symmetric n x n matrix, and the right-hand side is the com- 
mutator of matrices. With the help of N one can give Sym(n) the structure of a Lie algebra, 
the Lie algebra bracket being 



{wi,W2} = WiNw2 — W2NW1 



wi,W2 € Sym(n). 



Can we find a Lagrangian I £ C°°(Sym(n)) for which Equation (jlSh is of Euler-Poincare type 
with respect to the above Lie algebra? The answer is in fact given in [4j: a corresponding 
Lagrangian is 

ZH = itr(u;2). (19) 

We will show that the reduced Helmholtz conditions, applied to the current Lie algebra and 
dynamical system, lead to the correct Lagrangian. 

To make things more accessible we will consider only the case n = 2. For a basis of the Lie 
algebra we take the matrices 



Ex 



\ — ' 


■ 


, Ey - 


' 


1 " 








1 






and E^ 




1 



Further, without loss of generality we can take N to be 

1 
-1 

The non- vanishing Lie algebra brackets are then {Ex,Ey} = 2Ez, {E^, E^} = Ey and {Ey, E^} 
2Ez . An arbitrary element of the Lie algebra is of the form 

w = xEx + yEy + zE. 

and Equation (fT8|) is 



X y 
y z 



X 


y 




. y 


z 





-2y{x + z) 



X 



2y{x + z) 



We use now the notation of the Lie algebroid formulation of the Helmholtz conditions from 
Section 7. For $ we find 

^Ex) = {-3y^ + ^z^)Ex + {lxy-2yz)Ey + {4y^ -^xz)Ez, 
^{Ey) = {3xy-4yz)Ex + {4xz-^x'^ -^z^)Ey + {3yz-Axy)Ez, 
^{Ez) = {^y'^ -\xz)E^ + {lyz-2xy)Ey + {-3y^ + \x^)Ez, 



and for V 



VEx 

VEy 

VE, 



-{x + \z)Ey - yEz, 
{2x + z)Ex - {x + 2z)Ez, 
yEx + {z + \x)Ey. 
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The V-equations in this case (taking the symmetry of kij into account) are 



likxx) + (2a; + z)kxy + 2ykxz = 0, 

2z)kyy + ykyz " (2^ + z)kxx + (x — 2z)kxz — 0, 
li^xz) + (x + ^z)kyz + ykzz ~ ykxx ~ {z -\- \x)kxy = 0, 
j{kyy) - 2(2x + z)kxy + 2(x + 2z)A:j^^ = 0, 



li^yz] 
l{kzz) 

The ^-equations are 



(2x + z)kxz + {x + 2z)kzz - ykxy - {z + \x)kyy = 
2ykxz 



{2z + x)kyz = 0. 



kxxi^yx - Ayz) + kxy{-^x'^ + 4x2: - fz^) + kxz{-'kjx + Syz) 
-3y2 + iy2) _^ fcyj^dyx - 2yz) + kyzi^^y"^ - \xz), 



^xy\ 



kxxi^y"^ - \xz) + kxy{-2yx + fyz) + A;:r^(-3y2 + i^^) 

= kxzi-^y"^ + Iz"^) + kyzi^xy - 2yz) + A;^^(4y2 - ^xz), 
kxyi^y"^ - \xz) + kyy{-2xy + + kyziSy^ + ^x^) 



(20) 



(21) 



kxzi^xy - Ayz) + A:j^^(-|x2 



4xz - fz^) + /c2^(-4xy + 3yz). 



We will first try to find a solution of (j20p in which all are constants. In that case, adding 
one half times the (y, y)-equation to the (x,x)- and (z, z)-equations gives kxy = from which 
also kxz = and kyz = 0. Then the (x, z)-equation gives kxx = kzz and the (y, 2;)-equation gives 
kzz = \kyy. So the solutions of pop with constant coefficients are of the form 



1 
2 
1 



(22) 



It is easy to see that a multiplier k of this form also satisfies Equations (j2ip . The Hessian of the 
function 

l{x,y,z) = I(x2 + 2y2 + z2) 
takes the above form; this is exactly the Lagrangian (jl9p . 



An expression for the most general solution of Equations (|20p and (j2ip is beyond the scope of 
the current paper. However, instead of looking for constant solutions kij as above, we could use 
an additional symmetry assumption. For example, it seems natural to require that kxx — k^z 
(but that they are not necessarily constants). A tedious calculation reveals that in that case the 
only possible solution of the reduced Helmholtz conditions is again the multiplier in ()22p with 
constant coefficients. 



8.3 An illustrative example on the Lie group of the afRne line 

There are only two distinct Lie algebras of dimension 2. In this example we will use the Lie 
group of the affine line (the Euclidean group). An element of this group is an affine map 
R ^ R : t exp(qi)t + q2 and can be represented by the matrix 

exp(gi) q2 

1 ■ 
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The corresponding Lie algebra is given by the set of matrices of the form 

X y 

_ ■ 

A basis for this algebra is 





" 




' 


1 " 















for which {E^, Ey} = Ey. Let A = aE^ + bEy be a constant vector in the Lie algebra. We will 
determine whether there exists a regular Lagrangian for the dynamical system 

w = {w,{w,A}}, 

or, in the above basis, 

i; = 0, ij = x{hx — ay). 

For this system 

^{E^) = \{a- ifxyEy, VE^ = {-bx + i(a - l)y)Ey, 
^Ey) = -\{a- ifx^Ey, VEy = \{a + l)xEy. 

The V-equations are therefore 

x{hx - ay)kxxy - 2{-xb + ^{a - l)y)k^y = 0, 

x{bx - ay)kxyy -\{a + l)xkxy - {-bx + ^(a - l)y)kyy = 0, 

xibx — ay)kyyy — (a + l)xkyy = 0, 

and the only ^-equation is 

-(a - l)'^x'^kxy = (a - l)'^xykyy. 
If we differentiate the ^-equation with respect to x and y we obtain two more equations for the 

(o 1) (x kxxy ~\~ '^xkxy) — (o 1) {xykxyy + ykyy)i 

(ffl l)x'^ kxyy = (fl ^^'^ {-^y^yyy ~^ -^^yy)- 

The component kxx of the Hessian and its derivative kxxx are absent from these equations, and 
they will also not show up in any derived equation; there will therefore always remain freedom 
of choice for the x-derivative of kxx- We get 6 homogeneous linear equations in the 5 unknowns 
kxy, kyy, kxxy, kxyy and kyyy. If the rank of this system is less than 5 the system will have a 
non-zero solution. When a = 1, the rank is clearly 3. It can easily be verified that in all other 
cases the rank is 4. 

For reasons of clarity, we will deal first with the case where a = 1. 

1. The case where a = 1. In this case the ^>-equation is identically satisfied. The V-equations 
are now 

x(bx y)kxxy ~l" 'ibxkxy = 0, 
x(bx y)kxyy xkxy ~\~ bxkyy — 0, 
x{bx y^^yyy 2xkyy —— 0. 
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From the first and the last of these equations, we find that 

/i(y) 



k 



xy 



{bx - yf 



and 



k 



f2{x) 



yy 



{bx - yf 



respectively, as long as x ^ and bx — y 0. By substituting this result in the second equation 
and by interpreting kxyy once as dkyy/dx and once as dkxy/dy, we get a system of ODE's from 
which we can determine fi{y) and f2{x)- They are 



= -ba2-aiy and /2(x) 
The solution of the V-equations is therefore of the form 



aix + a2- 



k = 



2(6x-y)2 ^ J \X) 

ba2+aiy 

{hx-y)'^ 



ba2+aiy 
(bx—y)'^ 

aix+ag 
(bx-y)'^ 



det(A;) 



f{x){aix + 02) - al 
{bx - y)2 



This matrix is, however, not defined on the whole of R^. By continuity it exists on x = but 
it is not defined on bx — y = 0. So, there is no regular multiplier on R^. 

This is not the end of the story, however. The dynamical equations are now x = and y = 
x{bx — y). Notice that the lines x = and bx — y = ^ are both invariant under the flow. They 
divide the space R^ into regions, each invariant under the flow of the dynamical system. The 
matrix above is well-defined on the invariant region with bx — y ^ 0. It will be a multiplier 
provided its determinant is not zero, that is, provided f{x){aix + 02) — 7^ 0. The function 
l{x, y) = —a2 In \bx — y\ — cti-x In \bx — y| + a^y + h{x), with h"{x) = f{x), has the above matrix 
as its Hessian. However, for / to give the required Euler-Poincare equations, 02 and 03 must 
vanish. A non-degenerate Lagrangian on 6a; — y 7^ is therefore 

l{x, y) = — aixln \bx — y\ + h{x), 

where cti is a non-vanishing constant, and h{x) is an arbitrary function which is not of the form 
ai(xln \x\ — x) + a^x + as for any constants 04 and 05. 

In the following cases it will happen that there is no regular Lagrangian defined on the whole of 
R^, but it may be possible to find Lagrangians for subsets of R^ invariant under the dynamical 
flow. 

2. The case where a 7^ 1. In this case, the ^-equations come into play. As before, we can 
search first for the most general class of solutions of the V-equations, and then restrict to only 
those that also satisfy the ^-equation. Notice that e.g. the last of the V-equations leads to a 
further division of this case in subcases. We have 



kyy — < 



f2{x){bx - ay) 
/2(x)exp 
. 0, 

We will only summarize the results. 

2A. The case where a ^ 0. There is a regular Lagrangian of the form 



aj^O, 

a = 0, 6 7^ 0, 
a = 0, 6 = 0. 



l{x,y) = - — —\ay — bx\^ a\x\a+h{x), 
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where ai is a non-zero constant and h{x) is an arbitrary, but non-affine function. The Hnes 
X = and ay — bx = are invariant under the flow. 

2B. The case where a = and b 0. There is a regular Lagrangian of the form 



l{x, y) = aib^x exp j + h{x), 



where ai is a non-zero constant and h{x) is an arbitrary but non-affine function. The hne a; = 
is invariant under the flow. 

2C. The case where a = and 6 = 0. This is a degenerate case, there is no regular multiplier. 

Having decided in all cases whether a Lagrangian I on the Lie algebra g exists or not, it is 
instructive to give an expression for the corresponding Lagrangians L at the level of the Lie 
group G. If (q'i,(?2) are coordinates on the Lie group, then a left-invariant basis of vector fields 
is given by 

d - d 

Ex = Ey = exp(gi) — . 

oqi dq2 

Fibre coordinates (w*) = {x,y) with respect to this basis and (^1,32) with respect to the coor- 
dinate field basis are related as x = qi,y = exp (— (71)52. A right-invariant basis is 

Ex = -5 1- 92-^—5 Ey = ——. 

oqi Oq2 dq2 
The complete and vertical lifts of the left-invariant basis fields are 
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We can now rewrite a second-order field F in any of the following forms 

r = qi^ \-q2^ h/i^ + /2^ 

oqi dq2 dqi dq2 

= qi-R \- exp(-gi)g2 exp(gi) h qi-^ 



dqi \ \dq2 oq2 

d d 
+ fi-K^ + (exp(-gi)(/2 - qm)) exp(gi) — 
oqi oq2 

= xE^ + yE^+T.E^ + TyEl 
In the example under consideration, Fj, = and F^ = x{bx — ay), so 

/i = 0, /2 = (1 -a)gig2 + exp(gi)6g?). 

Let's look, for example, at case 2B (a = 0), where we have stated above that there exist a 
regular Lagrangian on the Lie algebra of the form l{x,y) = aib^xexp(y/bx) + h{x) {ai 7^ 0, /i 
non-affine) . By using left translations we can extend this to a Lagrangian on the whole of TG: 

T( ■ ■ \ u2- /exp(-gi)g2\ , , / • x 
L{qi,q2,qi,q2) = aib gi exp \^ — J +h{qi). 

Obviously, this Lagrangian is invariant: 

EnL) = ^ + q2^ + q2^ = and E^iL) = ^ = 0. 
oqi dq2 Oq2 Oq2 



28 



A short calculation shows that the Euler-Lagrange equations for the above Lagrangian do indeed 
return the differential equations qi = and q2 = Q1Q2 + beyip{qi)qf, as they should. 

Our analysis reveals only whether there is an invariant Lagrangian. In the case 2C (a = 6 = 0) 
where no such Lagrangian exists there could still be a (necessarily non-invariant) Lagrangian 
for the second-order system iji = 0, q2 = qiq2 on the two-dimensional Lie group. In [9] Douglas 
gave a more-or-less complete classification of the inverse problem for two-dimensional systems. 
A modern geometric approach to Douglas's classification can be found in [22]. A meticulous 
analysis using the methods described there shows that a regular Lagrangian must exist, even 
in the case 2C where we concluded that there is no invariant Lagrangian. In more detail, our 
case 1 belongs to Douglas's case I, and our cases 2A, 2B and 2C to his case Hal. 

Observe that if a = 6 = we are back in the example of the canonical connection. According to 
[23] . the most general Lagrangian for the case 2C, subject to the regularity condition, is given 
by 

L{qi,q2,qi,q2) = q^6{qi,q2,z) +tp{q^), z = 92/^1, 
where ip is an arbitrary function and is a solution of the PDE 

(subscripts denote derivatives, as usual). For example, the function 

L{qi,Q2,qi,q2) = ^qf + exp{-qi)-^ 

Zqi 

is a Lagrangian for the system in 2C. It is clearly not invariant since 

^f(L) = exp(-5i)|i. 

2qi 

In fact, there does not exist a function 9 for which the Lagrangian is invariant and regular. The 
relations Ei{L) = and ^2 (-^) = imply that 9q-^ + zOz = and Oq^ = 0, respectively. By 
taking the z-derivative of the first relation and by applying the second in the defining relation 
of 0, we can conclude that also 9z = ^ and 6*^^ =0. But then is a constant and the Lagrangian 
is clearly degenerate. 



9 Outlook 

We discuss briefiy two possible extensions of the current framework. First of all, let M be a 
manifold with a given symmetry group G. One can then set up an inverse problem for G- 
invariant Lagrangians on M. In that case, it has been shown in [5] that the Euler-Lagrange 
equations reduce to the so-called Lagran^e-Poincare equations. On the other hand, the technique 
of adapted frames can be easily extended to manifolds with a symmetry group; a description of 
the reduced equations for arbitrary second-order equations can be found in [8]. So the question 
would be when these reduced equations are of Lagrange-Poincare form. 

The second extension can be situated at the level of Lie algebroids. In this paper, we have 
discussed the inverse problem for Lagrangians on a Lie algebra Q. The original inverse problem 
deals with Lagrangians on TM. Both and TM are the two simplest cases of a Lie algebroid. 
So it seems natural to study an inverse problem for arbitrary Lie algebroids (the corresponding 
Lagrangian equations were given in e.g. [E]). The situation in the previous paragraph then 
coincides with the case that the Lie algebroid is TM/G, the so-called Atiyah algebroid. 
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